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1. Monopole and photon part of U(l) theory 

We are interested in the relationship between monopoles and classical chaos 



U(l) gauge fields are decomposed into a monopole and photon part across the phase 
transition from the confinement to the Coulomb phase. We analyze the leading 
Lyapunov exponents of such gauge field configurations on the lattice which are 
initialized by quantum Monte Carlo simulations. We observe that the monopole 
field carries the same Lyapunov exponent as the original U(l) field. As a long 
awaited result, we show that monopoles are created and annihilated in pairs as a 
function of real time in excess to a fixed average monopole number. 
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■ across the phase transition at (3 C ~ 1.01. We begin with a 4c? U(l) gauge 

theory described by the Euclidean action S{Ui} — (3 5^ p (l — cos9 p ), where 
Ui = U x ^ = exp^f?^) and 9 p = XjP + 9 X+ ^ V - 9 x+ p^ - 9 X , V . Follow- 
ing Ref. 1, we have factorized our gauge configurations into monopole and 
photon fields. The U(l) plaquette angles 9 X4iV are decomposed into the 
"physical" electromagnetic flux through the plaquette 9 X ^ V and a number 
nix,^ of Dirac strings through the plaquette 

where d x ,uv £ ( — T) + 71 "] an d t^x^v 7^ is called a Dirac plaquette. 

2. Classical chaotic dynamics from quantum Monte Carlo 
initial states 

Chaotic dynamics in general is characterized by the spectrum of Lyapunov 
exponents. These exponents, if they are positive, reflect an exponential 
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divergence of initially adjacent configurations. In case of symmetries inher- 
ent in the Hamiltonian of the system there are corresponding zero values of 
these exponents. Finally negative exponents belong to irrelevant directions 
in the phase space: perturbation components in these directions die out 
exponentially. Pure random gauge fields on the lattice show a characteris- 
tic Lyapunov spectrum consisting of one third of each kind of exponents. 2 
Assuming this general structure of the Lyapunov spectrum we present here 
its magnitude only, namely the maximal value of the Lyapunov exponent, 

The general definition of the Lyapunov exponent is based on a distance 
measure d(t) in phase space, 

L:= lim Urn -In M ( 2 ) 

t^ood(o)^o t d(0) 

In case of conservative dynamics the sum of all Lyapunov exponents is zero 
according to Liouville's theorem, "^2 Li — 0. We utilize the gauge invari- 
ant distance measure consisting of the local differences of energy densities 
between two 3c? field configurations on the lattice: 

d: =^E P \^Up-tvU' P \. (3) 

Here the symbol J2 P stands for the sum over all Np plaquettes, so this 
distance is bound in the interval (0,27V) for the group SU(N). Up and 
Up are the familiar plaquette variables, constructed from the basic link 
variables U x .i, 

U x ,i = exp (aA^T) , (4) 

located on lattice links pointing from the position x = (x\, x 2 , x 3 ) to x+aei. 
The generator of the group U(l) is T = —ig and A x .i is the vector potential. 
The elementary plaquette variable is constructed for a plaquette with a 
corner at x and lying in the ij-plane as U x ^j — U Xi iU x ^-ijU^. + j jU], j. It is 
related to the magnetic field strength B x y. 

U Xtij = exp (e ijk a 2 B x . k T) . (5) 

The electric field strength E x ,i is related to the canonically conjugate mo- 
mentum P Xi i — aU x .i via 

e ** = ^ (/ '''-'';.,) • ( 6 ) 
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The Hamiltonian of the lattice gauge field system can be casted into the 
form 



\(P,P) + l-\(U,V) 



(7) 



Here the scalar product stands for (A, B) — ^tr(AB^). The staple variable 
V is a sum of triple products of elementary link variables closing a plaquette 
with the chosen link U. This way the Hamiltonian is formally written as 
a sum over link contributions and V plays the role of the classical force 
acting on the link variable U. 

We prepare the initial field configurations from a standard four dimen- 
sional Euclidean Monte Carlo program on a 12 3 x 4 lattice varying the 
inverse gauge coupling (3. 3 We relate such four dimensional Euclidean lat- 
tice field configurations to Minkowskian momenta and fields for the three 
dimensional Hamiltonian simulation by selecting a fixed time slice of the 
four dimensional lattice. 



3. Chaos and confinement 




Figure 1. Exponentially diverging distance in real time of initially adjacent U(l) field 
configurations on a 12 3 lattice prepared at fi = 0.9 in the confinement (left) and at 
/3 = 1.1 in the Coulomb phase (right). 



We start the presentation of our results with a characteristic example of the 
time evolution of the distance between initially adjacent configurations. An 
initial state prepared by a standard four dimensional Monte Carlo simula- 
tion is evolved according to the classical Hamiltonian dynamics in real time. 
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Afterwards this initial state is rotated locally by group elements which are 
chosen randomly near to the unity. The time evolution of this slightly ro- 
tated configuration is then pursued and finally the distance between these 
two evolutions is calculated at the corresponding times. A typical exponen- 
tial rise of this distance followed by a saturation can be inspected in Fig. 1 
from an example of U(l) gauge theory in the confinement phase and in the 
Coulomb phase. While the saturation is an artifact of the compact distance 
measure of the lattice, the exponential rise (the linear rise of the logarithm) 
can be used for the determination of the leading Lyapunov exponent. The 
left plot exhibits that in the confinement phase the original field and its 
monopole part have similar Lyapunov exponents whereas the photon part 
has a smaller L max . The right plot in the Coulomb phase suggests that 
all slopes and consequently the Lyapunov exponents of all fields decrease 
substantially. 

4. Real-time evolution of monopoles 



P= 1-1 
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Figure 2. The number of monopole pairs as a function of real time evolution at (3 = 0.9 
(confined, left) and at /3 = 1.1 (Coulomb phase, right). 

An important subject of the present study is the evolution of monopoles in 
real time. 4 The outcome is displayed in Fig. 2 and shows that their number 
stays constant on average with creation and annihilation of a few monopole 
pairs only. The fluctuations decrease with increasing inverse gauge coupling 
P. A decrease of the total number of monopoles takes place towards the 
Coulomb phase. 

The real-time evolution of the monopole content in the decomposed 
U(l) field is depicted in Fig. 3. Although there are less monopoles in the 
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Figure 3. Monopole pairs in the decomposed field at /3 = 0.9 (left) and at j3 = 1.1 
(right). In the latter case the decomposition was repeated until complete absence of 
monopoles in the photon part. Full lines correspond to the monopole, dashed lines to 
the photon (residual) configuration. 



photon field than in the monopole field, their fluctuations are comparable 
in the confining phase. On the right plot, made for the Coulomb phase at 
(3 = 1.1, additionally the decomposition to photon and monopole fields has 
been repeatedly applied until the photon part had no monopole in it. This 
zero value is conserved during the real-time evolution. 
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